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ABSTRACT 

We consider the depletion of elements from the interstellar gas onto a 
population of very small dust grains. Adopting a grain model in which ~ 4% of 
the cosmic C abundance is in grains with radii < 10 A, we find that the rate 
of accretion onto these grains is adequately fast to account for the observed 
large depletions of elements like Ti, without invoking unreasonably high rates of 
mass transfer between interstellar phases or low grain destruction rates. If these 
grains are composed of arene rings, then only a limited number of metal atoms 
can be locked up in them. The depletion would be quenched when this limit is 
reached, unless there is a mechanism for transferring the metals to larger grains 
and refreshing the very small grain population, for example by grain coagulation 
and shattering in the diffuse ISM. If Fe depletes onto the very small grains, 
then for reasonable coagulation rates there is at least one metal atom per five 
C atoms in the very small grain population. Furthermore, ~ 60% of the cosmic 
Fe is associated with the carbonaceous grain population. It is unclear whether 
this scenario is compatible with observations. However, if there is another 
population of very small grains, with a large capacity for holding Fe atoms, it 
might be the sink for the most heavily depleted elements. 



Subject headings: dust — ISM: abundances 



1. Introduction 



Ultraviolet and optical absorption line studies have shown that many elements in the 
interstellar gas have lower abundances, relative to hydrogen, than in the Sun or in nearby 
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B stars. The depletion of species X is defined by 




gas phase 



cosmic 



(1) 



where n(X) is the number density of species X and riu is the number density of hydrogen. 
Savage & Sembach (1996) have recently reviewed depletion measurements for various 
interstellar environments. 

The conventional explanation for interstellar depletions is that some of the atoms 
have been accreted by dust grains. However, efforts to quantify this hypothesis for the 
most heavily depleted elements in diffuse gas have encountered difficulties. For example, 
Draine (1990) found that implausibly high rates of mass transfer from molecular clouds 
to the diffuse medium were required, but noted that the inclusion of a population of very 
small grains might help. The existence of such a population is inferred from observations 
of 3 to 60yum infrared emission, presumably generated by grains small enough to reach 
temperatures of 30 - 300 K or more upon the absorption of a single starlight photon (e.g. 
Draine & Anderson 1985). Very small grains are expected to be more effective accretors 
than bigger "classical" grains for two reasons. First, they contribute more total surface 
area. Second, they tend to be negatively charged or neutral, whereas the classical grains 
tend to be positively charged. Since the depleted metals in the gas are generally singly 
ionized, the collision cross sections for small grains are enhanced and those for big grains 
are diminished. 

In this study, we repeat the analysis presented by Draine (1990), but including a 
population of very small grains. In particular, we will aim to account for the depletion of 
Ti, which has been observed in both a cool diffuse cloud and a warm diffuse cloud along the 
line of sight towards ( Oph. For the cool cloud towards ( Oph, <5(Ti) = 1 x 10~ 3 ; for the 
warm cloud, S(Ti) = 5 x 10~ 2 (Savage & Sembach 1996). Although these particular clouds 
may not be typical, the severe Ti depletions appear to be. For sight lines along which 
the average H number density is 3 cm -3 , <5(Ti) = 2 x 10 -3 (Jenkins 1989), and Joseph & 
Jenkins (1991) found that 5(Ti) < 4 x 10 -4 for the cool clouds towards n Sco. Although Ca 
is more severely depleted than Ti towards ( Oph, Crinklaw, Federman, & Joseph (1994) 
conclude that generally Ca is depleted to the same extent as Ti, as a function of average H 
number density along the sight line. Depletion studies involving Ca are complicated by the 
fact that some Ca is expected to be doubly ionized, but Ca ++ has no resonance lines at 
wavelengths longward of 912 A. 

In §2 we derive expressions for the rate at which ions are removed from the gas by 
accreting onto grains, and compute accretion timescales for three interstellar phases. In §3 
we consider a scenario for mass exchange between phases and derive the constraints which 
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the observed Ti depletions impose on the timescales for dust destruction and mass transfer 
between phases. We conclude in §4 by discussing the plausibility of the resulting scenario. 



2. Depletion Rates 
2.1. Rate Coefficients 

Ions arrive at the surface of a spherical grain at a rate 

,oc /2\V2 / m x 3/2 / mv 2\ 

Rm = J dV U \m V eX P(^J^na>)™, (2) 

where k is Boltzmann's constant, T is the gas temperature, m is the ion mass, 6 max is the 
maximum impact parameter for which an ion strikes the grain, and n is the ion number 
density. When the grain is charged the ion moves in a Coulomb potential for which (Spitzer 
1941) 

( for raf < 

6max(u)=i 2^,^1/2 „ _, 2 z ^. e 2 , (3) 
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where Z g e is the grain charge, Z^e is the ion charge, m is the ion mass, and a is the grain 
radius. For charged grains, we neglect the modifications to 6 max which result from the 
"image charge" contribution to the interaction (Draine & Sutin 1987). Neutral grains are 
polarized by the ion's electric field. Treating the grain as a perfect conductor, Draine & 
Sutin (1987) found that 



b max (v) = a 
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mv 2 a 



1/2 



(4) 



The accretion timescale, r a , is defined and computed from 

-l 1 dn 1/2 (&ixkT\ l/2 f , ,(in firn/ 

T a 1 = = A 1/2 s / daa 2 — r-D(a), (5) 

ndt \ m p J J da y h w 

where A is the mass number of the ion, s is the sticking coefficient of the ion on the grain, 
trip is the proton mass, and n gr (a) is the number density of grains with radii less than a. 
The enhancement factor D(a) depends on Zj and Z g . With Zj fixed, 



D(a) = Y,f(Z 9 ,a)B(Z g ,a), 



(6) 
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where (Draine & Sutin 1987) 



B(Z g ,a) 



(l - Z ~§f) f ^ ^ < . (7) 
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The grain charge distribution function, f(Z g ,a), is obtained using the method described 
by Weingartner & Draine (1998). In Figure 1 we plot D(a), with Zi = 1, for graphite 
and silicate grains in three interstellar environments: cold neutral medium (CNM), warm 
neutral medium (WNM), and warm ionized medium (WIM). The adopted hydrogen number 
density n H , gas temperature T, and ratio of electron density to hydrogen density x e , for 
these three environments, and for molecular clouds (MC), are displayed in Table [I]. Since 
photoelectric emission is important in the grain charging, a radiation field must be specified; 
we used the spectrum of Mathis, Mezger, & Panagia (1983). The absorption cross sections 
were calculated using a Mie theory code derived from BHMIE (Bohren & Huffman 1983), 
with dielectric functions as described by Draine Sz Lee (1984) and Laor & Draine (1993). 
In molecular clouds, the grains are predominantly neutral and Zi = 0; thus D(a) ~ 1. 



2.2. Grain Population 

We consider two model grain size distributions, both of which are extensions of the 
MRN distribution (Mathis, Rumpl, & Nordsieck 1977), 

dn gr = Ca~ 3 ' 5 nnda, a min < a < a max , (8) 

with coefficients given by Draine & Lee (1984): C = 10~ 25 ' 13 cm 35 for graphite and 
C = 10~ 2511 cm 3,5 for silicate. The MRN distribution runs from a mm = 50 A to 
a max = 0.3/xm, but for the graphite component we extend the lower limit to 4 A, about the 
minimum size expected to be stable against sublimation in the interstellar radiation field 
(Guhathakurta & Draine 1989). We do not similarly extend the silicate component, because 
the observed IR emission from diffuse clouds does not exhibit the 10/im silicate feature 
(Mattila et al. 1996; Onaka et al. 1996). For our first distribution (i), we retain the MRN 
power law throughout the extended size range. If graphite grains have a density of 2gcm~ 3 , 
then the grains with a < 10 A (< 420 C atoms), contain ~ 2% of the cosmic C abundance. 
This is conservative, as attempts to model the observed IR emission (e.g. Desert, Boulanger, 
& Puget 1990) have suggested values two or three times bigger. Therefore, for our second 
distribution (ii), we adopt for the graphite component the form suggested by Draine & 
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Anderson (1985), which differs from (i) only for a < 30 A, where 

dn gr = C(30A)°V 4( Wa. (9) 

In this distribution, grains with a < 10 A contain ~ 4% of the cosmic C abundance. The 
quantity 

P( a ) = (L\ (10) 

V ; 1.48/ sdlna K ' 

indicates how grains of different size contribute to the overall accretion rate. In Figures 2 

and 3 we show P(a) for Z { = +1 ions and CNM, WNM, and WIM conditions. Clearly, the 

very small grains completely dominate the accretion. 



3. Consequences 

3.1. Model for Mass Exchange 

We consider mass exchange between the following three gaseous environments: warm, 
cold, and molecular. Hot ionized gas is neglected because of its small mass fraction. We 
also consider the removal of mass from the ISM during the star formation process, and 
the addition of mass to the ISM from evolved stars (stellar winds and supernovae) and 
from the infall of extragalactic material. For simplicity, we assume that the total mass 
and the metallicity of the ISM remain constant in time. To properly account for changes 
in metallicity with time, it would be necessary to consider the dependence of the dust 
abundance on metallicity and the distribution of newly formed metals over the interstellar 
phases. The fraction of mass removal due to star formation occurring in phase i is denoted 
Pi, where i = c, w, or m for cold, warm, or molecular, respectively; we assume p m = 1 and 
Pc = Pw = 0. We consider mass addition (both from evolved stars and infall) only into the 
warm and cold media, with the same relative contributions from the two sources in both 
media. Then the mass input can be treated as if it were derived from a single source, with 
Ti depletion <5 in ; if there is no Ti in the accreted extragalactic material, then <5 in has the 
same value as the depletion in the stellar outflow (relative to its enhanced metallicity). The 
timescale for injection of the entire mass of the ISM is denoted by r in , and the fraction of 
the input that goes into phase i is denoted by 

Ti is removed from the gas by accretion onto very small grains, with timescales 
r a> i. Since we wish to explain the most heavily depleted elements, we make the extreme 
assumption that atoms are only returned to the gas from grains when the grains are 
destroyed. We denote the timescale for the destruction of all the grains in phase i as r dji 
and the mass fraction in phase i as In the warm and cold phases, grains are destroyed 
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principally in supernova shock waves. The timescale for destroying all the grains in the ISM 
in these shocks is denoted t^', the fraction of this destruction occurring in the cold medium 
is denoted g c [f c T r ,l = and / w t^w = (1 ~~ 9c) T & 1 ]- We will suppose a small amount of 

additional grain destruction to take place in molecular clouds, with a timescale r^ m - 



3.2. Steady State Equations 



We suppose the following mass transfers occur: from warm to cold, cold to warm, cold 
to molecular, and molecular to warm, as indicated in Figure 4. We define the transfer 
timescales 7y as the time to transfer the entire mass of phase i to phase j. Assuming the 
depletions in each phase to remain constant yields 



1 "^+E^(4-* J ) + ^(*n-^) = ^, 



T d,j 



(11) 
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where Si denotes the Ti depletion in phase i (McKee 1989; Draine 1990). Constant mass in 
each phase yields 

We treat S c , S w , r ajC , r a , w , T^ m , r in , S in , q c , q w = 1 - q c , g c , g w = 1 - g c , T d>m , and the 
mass fractions /j as known quantities determined by observation (/j, S c , S w ), the adopted 
grain model (r ajC , r ajW , r ajm ), the model for mass injection (r; n , <5i n , g c , g w ), the relative 
rates for grain destruction in the different phases (g c , g w ), and the assumed rate for grain 
destruction in molecular clouds (r^ m ). Recall that we assume all star formation to occur in 
molecular clouds (p m = 1, p c = p w = 0). There are then five equations and six unknown 
quantities: the grain destruction timescale Td, the depletion in molecular clouds S m , and 
the four transfer timescales (r wc , r cm , r cw , r mw ). If Td is also specified, the solution for the 
remaining variables is as follows: 



S m = \-[(l-S w )+g c (S w -S c )} + ^ 

[Td T d ,m 



fv/S-w fcSc Si] 



fw (S w - S c ) 
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(13) 
(14) 

(15) 

(16) 
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and 



T~cw fc (<^w <^c) 



fc$c _ 9c (1 - 5 C ) _ <5 W - £ c ^ _ 1 - 5 m ^ ^ g c (<5 W - 5 in ) 



7d 



r d,: 



-i 



(17) 



The destruction timescale Td is constrained by requiring the above quantities to be 
positive. For example, we assume 5 W > 5 C ; then by demanding r wc > we find 



Td > g c (1 - 4) 



fc^c q c (<^in - 4) 



Requiring <5 m > implies 



, r/-| r \ , /e r M I /w^w /c^c /m <^ii 

Td < [(1 - <)w) + 9c {o w - 5 C )\ 1 



7"a,c ^d m 'ir 



(19) 



3.3. Exemplary Solutions 

The mass of the ISM is estimated to be ~ 5 x 1O 9 M , and the star formation rate (or 
rate of infall plus stellar outflow) is ~ 5M yr~ 1 (Timmes, Diehl, & Hartmann 1997 and 
references therein); hence we adopt r in « 10 9 yr. In Table |2| we show the adopted values 
of r ai (for each grain size distribution), and the accretion timescales were computed 
using (HD, with A = 48 and s = 1. For the warm and cold media we take Zi = +1, while for 
molecular clouds we take Zi = 0, as well as raising a min to 150 A. The adopted values for 
the mass fractions fi are rather uncertain, and not necessarily constant in time. 

In Table |3] we present three exemplary solutions. Case A uses distribution (i), with 
5[ n = 0.1, and cases B and C assume distribution (ii), with 5{ n = 0.1 (case B) and 1 (case 
C) . For each case the adopted value of Td lies at the center of the allowed range of values, 
which generally has a width <^ 10 7 yr. 

The results are quite insensitive to q c (the fraction of infall entering the cool phase), 
which we set to q c = 0.1, and to g c (the fraction of grain destruction taking place in the 
CNM), which we set to g c = 0.05. Only S m is sensitive to Td iin , which we set to 1 x 10 10 yr, 
and to r a m (although r cm and r mw begin to increase significantly when r a m gets so large 
that 5 m > 5 C , but this regime is unrealistic). For distribution («), = 1.4 x 10 9 yr when 
5i n = and increases with increasing <5 in ; above 5i n = 0.72, there is no solution (although 
if T in is increased to 1.4 x 10 10 yr, there are solutions up through 5 in = 1.0). For (ii), Td 
increases from 6.3 x 10 8 yr for 5j n = to 1.6 x 10 9 yr for 5- m = 1. 
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4. Discussion 

4.1. Overall Timescales and Depletions 

Jones and collaborators have made a series of models for grain destruction in shocks 
(Jones et al. 1994, 1996). They find that the smallest grains are the most resilient, since 
most destruction results from nonthermal sputtering off betatron-accelerated grains, and 
the larger grains are more effectively accelerated by the betatron mechanism. Less frequent, 
faster shocks efficiently destroy smaller grains through thermal sputtering. Furthermore, 
they find that shattering in collisions between larger grains produces large numbers of very 
small grains. Presumably these very small grains later coagulate to form bigger grains, 
maintaining the overall grain distribution. Since atoms which have accreted onto very small 
grains are incorporated into larger grains via coagulation, it would seem most appropriate 
for our purpose to adopt the timescale for destruction of the entire distribution. Jones et al. 
(1996) estimate the lifetime of carbonaceous grain material throughout all the phases of the 
interstellar medium, r d = 6 x 10 8 yr, which about equals the value adopted in §3 for case B, 
but is about a factor of 2 smaller than the values adopted for our other exemplary cases. 
The derived transfer timescales are in the expected range (McKee 1989). Thus, our case B 
example shows that accretion onto very small carbonaceous grains can occur fast enough to 
result in the depletions of even the most heavily depleted elements, if they stick perfectly. 

4.2. Sticking Coefficients 

The "sticking" coefficient s = P1P2, where p\ is the probability that the impinging ion 
(e.g. Ti + ) will be trapped on the grain surface long enough (^ 10 _12 s) to get rid of its 
kinetic energy and "thermalize" , and P2 is the probability that the atom will then remain 
on the grain for ^ 10 8 yr, until the grain is overrun by a supernova blastwave. This requires 
that the atom not be removed by either photodesorption (by a UV photon) or chemical 
reaction (with impinging H or O, in particular). For cold (T ^ 20 K) grains in cold gas 
(T ^ 100 K), it seems likely that p\ is of order unity for species other than the inert gases. 

The probability P2 is most likely highly variable from one species to another. We would 
require p 2 to be of order unity for Ti and other species that are strongly depleted, and 
p 2 <C 1 for species like N and S, which are not depleted. Photodesorption can keep a grain 
surface "clean" of species which are susceptible to photodesorption (Draine & Salpeter 
1979), since the lifetime against photodesorption by ambient starlight could be as short 
as ~ 3000 yr for an atom or molecule. Efficient photodesorption has been observed for 
UV-irradiated H 2 ice (Westley et al. 1995, 1996), but the photodesorption cross sections 
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for various metals at the various bonding sites on very small carbonaceous grains are not 
known. The small (carbonaceous) grains might consist of layers of fairly extended arene 
sheets. In this case, an atom located between sheets may have its bond(s) broken by a 
photon, but then transfer its excess energy to the grain through collisions, and promptly 
form new bonds. Thus, variations in P2 from element to element might arise from different 
cross sections for bond breaking or from different abilities to diffuse out of the grain once 
bonds are broken. 



Another consideration bearing on the sticking of metals regards the finite number of 
bonding sites in the very small carbonaceous grain population. Klotz et al. (1995) noted 
that in the most common PAH complexes, each metal is bonded by six p z electrons on the 
same ring, implying that rings adjacent to a metal-bonded ring are not available to bond 
other metal atoms. They concluded that there can be at most one metal per two rings, 
or 12 C atoms. However, there are many different possible bonding modes (Maslowsky 
1993), so we shall be more flexible in our criterion for what constitutes a "full" grain. 
Suppose fast accretion can be maintained so long as there are at most (3 metal atoms per 
C in the very small carbonaceous grains. (Klotz et al. would only consider (3 < 1/12.) 
Then, the very small grains become saturated when their metal content reaches f3 and 
must be replaced by a fresh population of very small carbonaceous grains which are largely 
metal-free, and therefore capable of accreting metals. We will refer to this as "refreshment" 
of the small grain population, and will let r re denote the timescale for this refreshment of 
the small carbonaceous grains. This refreshment would presumably occur via coagulation 
of very small grains onto classical grains, balanced by the shattering of classical grains to 
create "fresh" very small grains. If the coagulation takes place in molecular clouds and 
the shattering in the warm medium, then the refreshment rate t^ 1 < r^ 1 , the rate for 
exchange of mass from CNM to molecular clouds. In this scenario, rapid refreshment would 
require rapid mass exchange between the molecular clouds and other phases. Alternatively, 
both shattering and coagulation might occur in the CNM. In this case, we demand that 
Trc > Tcoag, where the timescale for grains to exit the very small grain population by 
coagulation is roughly approximated by 



4.3. Saturation of Very Small Grains with Metals 



coag 




kms 



v 



) 




3 



) 



(20) 



The evaluation is for grain distribution (ii). Interstellar turbulence and radiation pressure 
could conceivably lead to relative speeds between grains v ~ lkms -1 , although estimated 
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drift speeds resulting from anisotropic radiation tend to be about a factor of 50 smaller 
(Weingartner & Draine 1998). 

If there are 7 metal atoms per C atom in the entire carbonaceous grain population, 
then freshly shattered very small grains will presumably have ~ 7 metals per C as well. 
Thus, grains of size a remain unsaturated if 

rM<rr x (a) = ((3- 1 )^f\, (21) 

where iVc(a) is the number of C atoms in a grain and R s t{a) is the rate at which metals 
stick to the grain. We express r™ ax (a) in terms of its value for a = 10 A: 



D(a) (10 A 



E 7 FT" > ( 22 ) 



r--(a) D(10A) V a J ^ r™ ax (10 A) M 
where the sum is over all metals M which accrete onto the very small grains. For the CNM, 



C ax (10A) M = 11.6yr(/3- 7 ) 



^C^COS ( 



48X l/2 -i-i 



(23) 



J M 



where x cos is an element's cosmic abundance, relative to H. To approximate 7(M), the 
contribution to 7 from metal M, we assume that M is not at all incorporated into the 
carbonaceous grains during their formation processes. Balancing the addition of M 
(by accretion in the CNM, with timescale r a)C ) against the removal of M (during grain 
destruction, with timescale ra, and star formation, with timescale r m ), we find 

7 (M)=(£^) ^ (- + -)" (24) 

where m gra is the mass of the entire carbonaceous grain population, per H atom. 

In Table J| we display r™ ax (10 A) M /(/3 — 7), and 7(M) for several elements, adopting 



m 



5.3 x 10~ 27 g/H. The values of r a c to be used in equation ( p4f) are determined 
from equations (|18D and (|l^) by demanding that r d rj 6 x 10 8 yr, and are denoted r^ c q 
in Table 4; the values of (48/A) 1//2 s for use in equation (|23|) follow from these required 



values of r a , C ) assuming distribution (m). Some metals might accrete onto only the silicates, 
resulting in the growth of these grains, and not contribute to the need to refresh the 
very small carbonaceous grains. With (48/A) 1//2 s = 1, and adopting the graphite form of 
distribution (ii) for the silicates, except with lower cutoff at 15 A (10 A), r ac = 2 x 10 6 yr 
(8 x 10 5 yr). Thus Mg and Si could possibly only accrete onto the silicates. Iron is a 
marginal case; however, it would be strange if Fe didn't stick well to the very small grains 
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when Ti and Ni do. Excluding Mg and Si, but including Fe, we find 7 = 8.3 x 10 2 and 
T™ ax (10A)/(/5 - 7) = 2.9 x 10 8 yr- Setting r rc (a) > r coag ^ 1 x 10 7 yr, we find the minimum 
possible value for (3 for which grains with size a can be effectively refreshed. This ranges 
from w 0.35 for a = 4 A to 0.12 for a = 10 A. 

Thus, if Fe depletion is due to accretion onto very small carbonaceous grains, then 
Fe is a major component in the composition of these grains, and ps 60% of the cosmic Fe 
resides in the entire carbonaceous grain population. Although there appears to be no reason 
to exclude this possibility, it is certainly unexpected. The alternative, that Fe depletes by 
accreting onto silicate grains, is unattractive: (1) it isn't clear that the accretion rate onto 
the silicates is fast enough, and (2) one would expect Fe to stick to arene grains if Ti does. 

4.4. Another grain population? 

Perhaps there is a population of non-carbonaceous very small grains, onto which 
the most heavily depleted elements accrete. Tielens (1998) has recently inferred that the 
timescale for the return of Si from grains to gas is about an order of magnitude less than 
that for Fe, since Si is much less depleted in the warm phase. In computing our Table 4, we 
took the same for each species, resulting in different transfer timescales, as exemplified 
by r mw , the time for transfering gas from molecular clouds to the warm phase. Of course, 
r mw should be the same for each species. If we take r ajC = 3 x 10 6 yr for Si, appropriate 
for accretion onto an MRN silicate distribution extended down to 15 A, then we find 
Td = 1.4 x 10 s yr and r mw = 7 x 10 7 yr. The same r mw is obtained for Ti and Fe if their 
accretion timescales are increased by slightly less than a factor of two, and = 1.0 x 10 9 yr. 
Thus, we can imagine a silicate grain population deficient in Fe and with a destruction 
timescale somewhat less than the 4 x 10 s yr derived by Jones et al. (1996), along with a 
more resilient, Fe-rich grain population which includes very small grains. Metallic Fe grains 
are ruled out because they would produce excessive magnetic dipole emission at 90 GHz 
(Draine & Lazarian 1999). Sembach & Savage (1996) have suggested oxides as a grain 
component, on the basis of gas phase abundance measurements in diffuse halo clouds. They 
found more Fe plus Mg per Si in dust than expected for pure olivines or pure pyroxenes. 

A 20/iin emission feature has recently been detected in emission in warm circumstellar 
dust shells, and attributed to FeO grains (Waters 1998). Since single-photon heating of 
the very small interstellar grains results in emission at A ^ 5/im, candidate materials 
will be constrained by the observed emission spectrum of diffuse clouds. Koike et al. 
(1981) and Henning et al. (1995) have experimentally determined IR optical constants for 
well-known iron oxides, and reported Q^bs/a for spheres in the Rayleigh limit (Q a bs is the 
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absorption efficiency factor). We reproduce their band positions in Table [5[ It is important 
to note that the band positions depend on the particle shape. Henning et al. (1995) also 
calculated absorption for two continuous distributions of ellipsoids and found that the FeO 
band position shifts from 19.9/xm to 21.0/xm or 23.4/zm. Serna, Ocana, & Iglesias (1987) 
and Wang, Muramatsu, & Sugimoto (1998) discussed the effect of particle shape on the 
hematite spectrum. We suspect that the 9/im feature in the Koike et al. hematite spectrum 
is due to an impurity in their sample; this feature is absent in the Nyquist & Kagel (1971) 
spectrum and present, but attributed to an impurity, in the Sadtler spectrum (Ferraro 
1982). Thus, the Matilla et al. (1996) and Onaka et al. (1996) studies do not rule out 
very small grains with iron oxide composition, since the longest wavelength they observed 
is 11.7/im. Spectrophotometry of the emission from diffuse clouds in the 15 to 30yum range 
would resolve the issue, but are beyond current capabilities. Thus at this time we cannot 
exclude the possibility that a significant fraction of interstellar Fe resides in ultrasmall Fe 
oxide particles. 

4.5. Conclusion 

Observations of infrared emission in the 3 to 60/zm range indicate that there is a 
population of very small dust grains, extending down to sizes of a few A. The presence of 
several spectral features associated with polycyclic aromatic hydrocarbons and the absence 
of the 10/im silicate feature imply that there are carbonaceous grains with a < 10 A, but 
no such silicate grains. The rate at which gas cations collide with these ultrasmall grains 
greatly exceeds that for collisions with bigger "classical" (a ^ 0.01/zm) grains. If the 
elements which are most heavily depleted in the ISM stick efficiently upon collision, then 
accretion onto very small grains occurs quickly enough to account for their depletions. 
This fast accretion is completely independent of the model adopted for the classical grains, 
which is currently controversial (Mathis 1996). Efficient sticking requires resistance against 
photodesorption in the interstellar radiation field and against chemical attack. Furthermore, 
if the depletion of the more abundant element, Fe, is due to accretion onto very small 
carbonaceous grains, then Fe would be a major component in the composition of these 
grains, at a level which might be chemically impossible. Perhaps Fe doesn't stick well to 
the very small carbonaceous grains, instead depleting onto a silicate population extending 
down only to ~ 10 to 15 A. Alternatively, there might be another population of very small 
iron-rich grains. There are well-known oxides which cannot be ruled out to date, but which 
have spectral features in the 15 to 30/mi range. 
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This preprint was prepared with the AAS IAT^X macros v4.0. 



Table 1. Idealized phases for interstellar matter 



quantity CNM WNM WIM MC 



%(cm- 3 ) 30 0.4 0.1 1000 

T(K) 100 6000 8000 20 

x e 0.0015 0.1 0.99 
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lO" 7 10" 6 10- 5 

a (cm) 



Fig. 1. — The enhancement factor, D(a), for the collision of ions with — +1 with grains. 
For graphite grains in the CNM, D drops to 4 x KT 12 at a = 3 x 10 cm. 
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lO" 7 10" 6 10- 5 
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Fig. 2. — -P(a), showing the contribution to depletion as a function of grain size, for the 
CNM. The discontinuity at 50 A is due to the population of silicate grains with a > 50 A. 
The lower curve is for grain size distribution (i) and the upper curve for distribution (ii). 
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Fig. 3.— Same as Figure 2, for the WNM and WIM. 
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Fig. 4. — Mass transfers between warm, cold, and molecular phases. Each arrow is labelled 
by the mass transfer rate divided by the total mass of the ISM. 
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Table 2. Input Values 



quantity 


w 


c 


m 


Vi (y r ) (l) 


1.1 x 10 8 


4.4 x 10 5 




r a ,i (yr) (a) 


6.0 x 10 7 


2.1 x 10 5 




r aji (yr) (MRN) 






5.6 x 10 6 


/i 


0.2 


0.3 


0.5 


^i(Ti) 


5 x 10~ 2 


1 x 10~ 3 





Note. - - For grain size distribution (i), the adopted value of r ajW is intermediate between 
the values of 1.0 x 10 8 yr for the WNM and 2.8 x 10 8 yr for the WIM. For distribution (ii), 
the corresponding values are 4.9 x 10 7 yr for the WNM and 1.5 x 10 8 yr for the WIM. 



Table 3. Exemplary Solutions 



case 


dist 


5m 


Td 




7~cm 


Tew 




5ra 


A 


(i) 


0.1 


1.4 x 10 9 


1.5 x 10 7 


4.5 x 10 7 


4.7 x 10 7 


8.7 x 10 7 


5.9 x 10~ 4 


B 


(ii) 


0.1 


6.3 x 10 s 


7.3 x 10 6 


2.4 x 10 7 


2.0 x 10 7 


4.3 x 10 7 


6.1 x 10" 4 


C 


(ii) 


1.0 


1.6 x 10 9 


7.6 x 10 6 


2.4 x 10 7 


2.1 x 10 7 


4.4 x 10 7 


6.1 x 10" 4 



Note. - For 5 W = 5 x 10 2 , 5 C = 1 x 10 3 , sticking coefficient s — 1, mass number A 
= 48 (Ti), q c = 0.1, g c = 0.05, r in = 1 x 10 9 yr, and r djm = 1 x 10 10 yr. All timescales are in 
yr. 



- 21 - 



Table 4. Required Refreshment Timescales 



M 


log x cos 










r rcq 

a,c 


r™ ax (10A) 
P-1 


7 (M) 






Mg 


-4.42 


0.13 




2.8 x 10" 


-2 


6.0 x 10 6 


3.1 x 10 s 


7.6 x 10- 


-2 


9.0 x 10 7 


Si 


-4.45 


0.30 




4.9 x 10- 


-2 


1.3 x 10 7 


4.1 x 10 s 


5.6 x 10- 


-2 


3.3 x 10 s 


Fe 


-4.49 


5.6 x 10" 


-2 


5.4 x 10" 


-3 


1.0 x 10 6 


3.2 x 10 s 


7.4 x 10- 


-2 


4.4 x 10 7 


Ca 


-5.66 






1.9 x 10- 


-4 


5.5 x 10 4 


7.3 x 10 9 


3.3 x 10" 


-3 




Ni 


-5.75 


3.1 x 10" 


-2 


1.8 x 10- 


-3 


3.4 x 10 5 


5.9 x 10 9 


4.1 x 10- 


-3 


2.9 x 10 7 


Cr 


-6.32 


8.5 x 10" 


-2 


5.2 x 10- 


-3 


1.0 x 10 6 


2.2 x 10 10 


1.1 x 10- 


-3 


6.6 x 10 7 


Ti 


-7.07 


5.0 x 10- 


-2 


1.0 x 10- 


-3 


2.0 x 10 5 


1.3 x 10 11 


1.8 x 10" 


-4 


3.7 x 10 7 



Note. The second, third, and fourth columns are taken from Savage and Sembach 
(1996), Table 5. All timescales are in yr. The depletions 5 C and 5 W are for the line of 
sight to ( Oph. r™ c q is the accretion timescale required to yield the observed depletions for a 
destruction timescale Td = 6 x 10 8 yr. The depletion 5 W for Ca in the warm cloud is unknown. 
We assumed that (1 — 5 W ) /(l — 5 C ) is the same for Ca as for Ti in computing r^ e c q . r™ ax (10 A) 
is the maximum allowed refreshment timescale for a = 10 A carbonaceous grains, where (3 
is the maximum allowed number of metal atoms per carbon atom, and 7 is the number of 
metal atoms per carbon in "fresh" carbonaceous grains. 



Table 5. IR bands for Fe oxides 



mineral 


IR band position 


reference 




/iiii 




hematite (a-Fe 2 03) 


9, 18, 21, 30 


Koike et al. 1981 


magnetite (Fe 3 04) 


17, 25 


Koike et al. 1981 


wiistite (FeO) 


19.9 


Henning et al. 1995 



Note. - - We suspect that the 9/zm feature in the hematite spectrum is due to an impurity. 



